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, Andrews-Baxter- Forrester models in regimes III and IV are calculated with fixed 

boundary conditions. The interfacial tensions are calculated between arbitrary 
| phases and are shown to be additive. The associated critical exponents are given 

by 2 — a s = u = v with v = (L + l)/4 in regime III and 4 — 2a s = u = v with 
v = {L + l)/2 in regime IV. Our results are obtained using general commuting 
transfer matrix and inversion relation methods that may be applied to other 
solvable lattice models. 
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1 Introduction 

There has been a recent convergence of interest in statistical mechanics on systems 
with a boundary. It is well known that two-dimensional lattice models without 
a boundary are exactly solvable [BJ by commuting transfer matrix methods if the 
local statistical weights satisfy the celebrated Yang-Baxter equation. It has also 
been known for some time, from the work of Cherednik g] and Sklyanin Q, that 
this integrability extends to vertex models defined on a strip with open boundaries 
provided the local boundary weights satisfy an additional reflection equation or 
boundary Yang-Baxter equation. More recently, reflection equations for interaction- 
round-a-face (IRF) models have been introduced H, ||, [|, and integrability has been 
established for lattice spin models defined on a strip with fixed or more general 
boundary conditions. 
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Once integrability with a boundary has been established there are various quan- 
tities of physical interest, such as the surface free energies and interfacial tensions, 
that one would like to calculate, and methods to achieve this need to be developed. 
In this direction surface free energies of the Andrews-Baxter-Forrester (ABF) mod- 
els in regime III |7|, ||], the eight- vertex model [[|, the dilute Al models |lO|| , and 
the CSOS models [11] have been obtained by an extension of the inversion relation 



method [12, [D| used to calculate the bulk free energies. 

In this paper we extend the analysis of J?]] to obtain the surface free energy of 
the ABF models in regime IV. In addition, we extend the generalized inversion 
relation method to calculate also the interfacial tensions and correlation lengths. In 
this way, we establish that for solvable lattice models it is possible to obtain the 
bulk free energies, the surface free energies, the interfacial tensions, the correlation 
lengths and their associated critical behaviours all by studying the relatively simple 
inversion relations. 

The layout of the paper is as follows. For the rest of this section we follow 
reference in describing the ABF models with fixed boundary conditions. In 
sections 2 and 3, we obtain the bulk and surface free energies. The band structure 
and ground state degeneracies of the eigenvalue spectra of the transfer matrices are 
discussed in section 4. The calculations of the interfacial tensions and correlation 
lengths are then given in sections 5 and 6 respectively. We conclude with a brief 
discussion. 

1.1 ABF models with fixed boundaries 



The ABF models [14| are restricted solid-on-solid models in which heights on the sites 
of the square lattice take values in the set {1, 2,3, ... ,L} subject to the condition 
that the values of heights on adjacent sites must differ by ±1. The Boltzmann 
weights depend on a crossing parameter A = vr/(L + 1), and a spectral parameter u. 
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The i?i(u) = $i(u,p) is a standard elliptic theta function with nome p. The tem- 
perature variable t = p 2 measures the deviation from criticality. The critical limit 
of the ABF models is t — > 0, approached from t > in regime III and t < in 
regime IV. We therefore express the nome p in terms of a real parameter e > as 

{ e -7rE regime III 
(1.4) 
ie n£ regime IV 

so that t = ±exp(— 2ire), and t —* as e —* oo. The product expansions of the 
functions t?i and #4 are given by 



t?i(u,p) = 2p x / 4 sinn J](l-2p 2n cos 2u+p 4n )(l- p 2n ) (1.5) 
n=l 

00 

4 (u,p) = JJ(l-2p 2n - 1 cos2n + p 2 ( 2n - 1 ))(l-p 2n ). (1.6) 



n=l 



The i?i functions satisfy the quasiperiodicity properties 



t?l(u + 7T,p)=-1?l(tX J p) (1.7) 

0i(t* -ihxp,p) = -p- l e- 2iu $i{u,p) (1.8) 



and the "conjugate modulus" transformations 



01 (u, e~™) = e -(«-V2) 2 /^ E ( e -2u/e je -2*/e) ( L g) 

0l(«, ie-^) = e W8 e -(«+V4) 2 M jB ( e «/^ _ e -V2^ ( L10 ) 

v 2e 

00 

£(a;,p) = JJ(1 -p™- 1 ^! -p^aT 1 )^ -p n ). (1.11) 



n=l 



Following reference ||] we introduce boundary weights. These weights depend 
on an additional real parameter £, which is independent of u and may be different 
for the left and right boundaries. The non-zero boundary weights are 

K ( o\ ^ ^!(( a ±l)A) V/ 2 ^(ni^^ATQ 

From the face weights and boundary weights we construct a double-row transfer 
matrix D(u). For a lattice of width N, the entry of the transfer matrix correspond- 
ing to the rows of heights a = {a\, . . . , qat+i} and b = {bi, . . . , bjv+i} is defined 
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diagrammatically by 



fei 61 b2 63 6jv b?f+i 6jv+i 




The solid heights {ci,... , c^+i} are summed over. As the boundary weights are 
diagonal, we must have a\ = b\ and a^ + \ = 6jv+l- Furthermore, these boundary 
heights, which we will call a L and a R , are fixed to the same values for all entries in 
the transfer matrix. The parameters £ L and £ R are similarly fixed for all entries. 
Defined in this way, the double-row transfer matrix exhibits the crossing symmetry 

D(\-u) = D(u) (1.13) 

and consequently is real symmetric for u real. More importantly, however, the 
double-row transfer matrices form a commuting family, 

D(u)D{v) = D{v)D(u). (1.14) 

This implies that the eigenvectors of D{u) are independent of u, so that functional 
equations satisfied by the transfer matrix are also satisfied by its eigenvalues. In 
particular, all eigenvalues of the transfer matrix satisfy the crossing symmetry (1.13). 



It should be emphasized that all the matrices in a commuting family share the same 
boundary heights a L and a R , and the same values of £ L and £ R . 

To ensure that the largest eigenvalue of the double-row transfer matrix is non- 
degenerate for all < u < A, we impose the restriction 

^<|e L |,|e R |<A (1.15) 

and in addition require that £; L,R > when a L ' R = 1 and £ L,R < when a L,R = L 
(note that the £ L ' R > and £ L,R < regions are disconnected). Proof of the 
sufficiency of these restrictions proceeds as follows. 

We first show that all the off-diagonal elements of D(u) are non- negative. Con- 
sider an off-diagonal element (a\D(u)\b) , with a and b such that (a\D(u)\b) is 
not identically zero for all u. Since the element is off-diagonal, we must have 
Oj 7^ hi for some 2 < i < N . We use the boundary crossing relation (equation 
(3.8) of reference Q) to replace the left boundary weight K{\ — u) with K(u). 
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This introduces a face with spectral parameter 2u — A and also a (positive) factor 
of $i(A)/i?i(2u). We then use the Yang-Baxter equation (equation (3.4) of refer- 
ence J3J) to push this face to the right until it separates the (i — l)-th and i-th pairs 
of faces. The upper and lower heights bi and a» differ, so from equation ([O]) the 
weight of this face is $i(2A — 2u)/$i(A), which is positive for < u < A. Since the 
other face weights are positive, and since the boundary weights K(u) are positive 
for < u < min {|£ L | , |£ R |}) the entry (a\D(u)\b) is the sum of positive terms. The 
crossing symmetry ( |1.13j ) then dictates that what holds for u must hold for A— u, and 
since A/2 < min {|£ L | , the off-diagonal elements of D{u) are non- negative for 

all < u < A. 

We now observe that the elements of D(u) are bounded, so there exists a real 
number M > such that all elements of the matrix MI + D(u) are non-negative for 
< u < A. The largest eigenvalue of this matrix is non-degenerate by the Perron- 
Frobenius theorem, and it follows immediately that the largest eigenvalue of D(u) 
is also non-degenerate. 

Indeed when u satisfies \u — A/2| < min { |£ L | , |£ R | } — A/2, the double-row trans- 
fer matrix is non-negative definite. This is due to the fact that when each of the 
boundary weights K(u) and K(X — u) is non-negative, D{u) is expressible as the 
sum of non-negative definite matrices. 

When u = A/2, the symmetry of the face weights is such that the model is 
isotropic. In this case the values £ = ±A/2 deserve special mention since for these 
choices the isotropic model has all boundary heights fixed. This is easily seen from 
the definition ( |1.12j ) as, for fixed a, only one of the choices a ± 1 gives a non-zero 
boundary weight. The non-zero boundary weights then contribute only a constant 
factor to each entry of the transfer matrix. Aside from this trivial factor, the lattice 
exhibits pure fixed boundary conditions, with boundary heights alternating either 
{a, a + l,a,a + 1, . . .} or {a, a — 1, a, a — 1, ... }. 

If we divide each of the face weights by $i(u — A/2), and each of the boundary 
weights by — A/2) 2 , then the quasiperiodicity ( |Q| ) implies that replacing u by 
u — i In p simply introduces a gauge factor to each of the normalized weights. These 
gauge factors cancel in the entries of D(u), so it follows that the quasiperiodicity of 
D{u) is that of $i(u — \/2) 2N+A . If we therefore define a normalized transfer matrix 



^(u-A/2) 2 



D{u) = , Q ,., I', (1-16) 
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then the entries and eigenvalues of D{u) are doubly periodic functions of u with 

{e regime III 
(1.17) 
2e regime IV. 

In regime IV there is an additional symmetry within the period rectangle 

b{u±Ti/2 + me) = D{u). (1.18) 

As was shown in [0], the eigenvalues of the ABF models with fixed boundary 
conditions satisfy the inversion identity 

S- lSl D(u)D{u + A) = e L e R s_ 2S2 /_i/i + slfoD 1 ' 2 ^) (1.19) 

where the functions and fk are given by 



M^ + k\) N 



n 27V 



(1.20) 



^i(A) 

The functions e L and e R depend on the left and right boundary conditions and, for 
generic a and £, take the form 



The function D 1,2 (u) is an eigenvalue of the double-row transfer matrix at fusion 
level 1x2. We observe as in the periodic case that for large N, foD 1,2 is 
exponentially small in —A/2 < Re(u) < A/2. Hence, for N large, the second term 
in the inversion identity (|1.19|) can be neglected. So in calculating bulk and surface 
properties we just need to solve the simple inversion relation 

s- lSl D{u)D{u + A) = e L e R s„ 2 s 2 /_i/i (1.22) 

and the crossing relation 

D{u) = D{\-u) (1.23) 

subject to the appropriate analyticity and quasiperiodicity in an open strip contain- 
ing < Re(u) < A. These properties determine the bulk and surface quantities 
uniquely. 

The largest eigenvalue D max (u) factorizes into contributions from the bulk, the 
surfaces and the interface 

) ~ Kb(u) 2N n s (u) A(u) as N — > oo. (1-24) 
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When the left and right boundary conditions favour the same phase, which is the 
case when n(a, (5) = in the notation of section 4, there is no interface and A(u) = 1. 
When the left and right boundary conditions favour different phases, which is the 
case when n(a,(3) > 1, there is an interface and the factorization applies to all the 
eigenvalues in the first band. Clearly, the inversion and crossing relations factorize 
into bulk and surface terms and can be solved sequentially for the bulk free energies, 
the surface free energies (n(a,/3) = 0) and the interfacial tensions (n(a,/3) > 1). 
This we do in the following sections. 



K h (u)n h {u + A) = _ a , M2 • (2.1) 



2 Bulk free energies 

Equating the bulk terms in the inversion relation ( |1.22j ) implies that the bulk par- 
tition function per site satisfies the functional equation 

fli(A -u)fli( A + «) 

This equation is of course independent of the boundary conditions on the lattice 
and agrees with the equation for periodic boundaries. It has been solved previously, 
but we include the solution here for completeness. We use the standard techniques 
developed by Stroganov |[^] and Baxter Ji~3| ], 

Since the eigenvectors of D(u) are independent of u, it follows that the eigen- 
values D(u) have the same analyticity and quasiperiodicity as the elements of the 
transfer matrix. In particular, the eigenvalues must be entire functions of the spec- 
tral parameter u. They are therefore completely characterized by their zeros and 
growth at infinity. In regime III the zeros of the largest eigenvalue accumulate on 
the lines Re(u) = —A/2 and Re(u) = 3A/2. In addition to these lines, in regime IV 
zeros accumulate on Re(u) = (3A — vr)/2 and Re(u) = (ir — A)/2 in accordance with 



the periodicity ( 1.18| ). The strip —A/2 < Re(u) < 3A/2 is free of order N zeros in 



both regimes. This is confirmed by numerical studies for large finite N. So inside 
this strip ^(u) is non-zero and ln«b(^) is analytic. The quasiperiodicity 

{ e j(A-2u) e 7r£* Kh ( u \ regime III 

V ' B (2.2) 
e 2i(\-2u) e 2ne* Kh ( u j regime Ty 

implies that the second derivative of lnKb(^) is periodic, with period ire*. Hence 
lnKb(tt) can be expanded in the form of a generalized Fourier series 

oo 

]n Kh (u) =Au 2 + Bu+ c k e 2ku/£ \ (2.3) 
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To evaluate the coefficients Ck, A and B, we take the logarithm in ( |2.1[ ), expand 
the right hand side using the appropriate conjugate modulus transformations (|L 
and Ql-lCp , and equate coefficients. 

2.1 Regime III 



In regime III we have t > 0, so we rewrite the right hand side of (2A) using the 



conjugate modulus transformation (|1.9[). With both sides of (|2.1[) expanded in 



powers of exp(2n/e), we match coefficients and impose the crossing symmetry ( 1.23 ) 
to obtain the solution 

1 cosh[(7r — 2X)k/e] sinh[(A — u)k/e] smh.(uk/e) 
lnKb(u) = " U)U + 2 ^ fcsinh(vrfe/ £ )cosh(Afc/ £ ) ' 

(2.4) 

Inside the region —A/2 < Re(u) < 3A/2 this function gives the bulk behaviour of the 
largest eigenvalue of the transfer matrix. Applying the Poisson summation formula 
to the infinite sum gives the behaviour of the free energy in the critical limit t —* + . 
When L is even lnKb(w) is regular, but when L is odd the leading-order singularity 
is 







In K b ~ t* /2X \nt. (2.5) 
Since In Kb ~ t 2 ~ a , when L is odd the bulk critical exponent a is given by 

2-a = 4 = ^±±. (2.6) 
2A 2 y J 

2.2 Regime IV 

In regime IV the temperature variable t is negative, so we use the conjugate modulus 
transformation ( |1.10| ) to rewrite the right hand side of ( |2.lD in powers of exp(n/e). 
As before, we match coefficients and use the crossing symmetry ( |1.23j ) to obtain 

, . 1 . cosh[(7r/2 — 2A)Ar/el sinh[(A — u)k/e] sinhfuk/e) 

in« b (it) = — (\-u)u+y — — ; — ' \ , , — w /, : — — 

K ' ire y 1 ^ A:sinh(7r/c/2e)cosh(Afc/e) 

Esinh[(7r/2 — 2X)k/e] sinh[(A — u)k/e] sinh(n/c/e) 
/tcosh(vrA:/2e) cosh(Afc/e) ' ^ ' ' 

fc+|=l 

Once again we apply the Poisson summation formula to obtain the leading order 
singularity [lj] 

-tW 2A L even 

In K b ~ { " (2.8) 
-t)^/ 2A ln(-t) L odd 
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so that in regime IV, the critical exponent a is given by equation ( |2.6| ) for all L. 



3 Surface free energies 

In order to calculate the surface free energies, we assume that the boundary condi- 
tions on the left and the right edges of the lattice favour the same phase, so that 
n(a, 0) = in the notation of section 4. In this case there is a unique largest eigen- 
value separated from the other eigenvalues by a gap. We can then divide out the 



bulk quantities in the inversion relation ( 1.22 ) to obtain an inversion relation for the 



surface partition function per site k s (u). Explicitly, 

, . . 4S ??i(2A-2u)??i(2A + 2u) T . . n . . ,„„, 

The form of this inversion relation suggests a natural factorization of k s (u) into a 
term independent of the boundary conditions, a term dependent on the left boundary 
condition, and a term dependent on the right boundary condition. We therefore 
write 

«b(«) = «° («)/<£ (u)/ef(t*). (3.2) 



The solution of equation (3.1) proceeds in a similar fashion to the solution of the bulk 
inversion relation, but whereas the analyticity of In Kb(u) depended on the absence 
of order TV zeros, the analyticity of ln« s (ii) depends on the absence of order 1 zeros. 
Our numerical studies show that when n(a, (3) = the largest eigenvalue of the 
double-row transfer matrix is indeed non-zero for < Re(it) < A. We therefore 
conclude that \jxk s {u) is analytic on this strip, and along with the quasiperiodicity 

e 4i(A-2n) e 47re* ^ ^ re gime III 



k s (it + ine*) 

this allows us to expand hiK s (u) as a generalized Fourier series. 



(3.3) 

e 8l ( A " 2n ) e 8ne * k s (u) regime IV 



3.1 Regime III 

In regime III we use once again the conjugate modulus transformation (|1.9|) and the 



crossing symmetry (1.23) to match Fourier coefficients and obtain 
ln^(n) = 

A sinh[(7r — 3X)k/e] sinh(AA;/e) cosh[2(A — 2it)A;/e] , . 

— C'T- )+ 2^ fcsinh(vrA:/e)cosh(2AA : /e) ( ' ' 

k=l 
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and for generic a and £, 
1 



ln^(u) 

oo 

fc=i 



7Tc 



[(aA + 0(t - 2 + " 2A)vr + ( 2 " a 2 )* 2 + 2 ^(A - u)] 



cosh[(aA + | - |£|)fc/e] cosh[(7r - aA - £ - |£|)fc/e] cosh[(A - 2tt)fc/e] 

sinh(-7rA;/e) cosh(Afe/e) 

+2 f; co f[ ( l7_ 2 r A y ] . (3.5) 



k=l 



fcsinh(7r/c/e) 



We note that, as it should, the height reversal transformation a —* L + 1 — a and 
£ — ► — £ leaves ( jO|) unchanged. 

When the Poisson summation formula is applied to the above expressions for 
hiKg(u) and lnK^' R (ii), the dominant behaviour as t — > + is seen to come from 
equation ( |3,4| ). If L = 2 (mod 4), however, this term is regular and we must consider 
the sum in equation fl3.5| ) . This latter sum is also regular when £ > with a odd and 
when £ < with a even Q, but such exceptions aside, the leading-order singularities 
of the surface free energy have the form 



hiK s ~ < 



^/4A L = or 1 (mod 4) 
t ir/2X L = 2 (mod 4) 
W 4A lni L = 3 (mod 4). 



(3.6) 



Since ln/t s ~ i 2 Qs , the surface critical exponent a s [17, 111] in regime III is given by 



(L + l)/2 L = 2 (mod 4) 
(L + l)/4 otherwise. 



(3.7) 



3.2 Regime IV 



In regime IV we use the conjugate modulus transformation ( l.lOj ) to match coeffi- 
cients and obtain 

, 0/ , A ^ sinh[(vr/2 - 3A)/c/e] sinh(A£;Ar) cosh[2(A - 2u)k/e] 
In K s {u) = — (tt-6A) + 2^ 



2TT£ 



k=l 



ksmh(irk/2e) cosh(2AA;/e) 



^ cosh[(vr/2 - 3A)fc/e] sinh(Afc/e) cosh[2(A - 2u)k/e] 
+ 2^ f.„„„u/_f./o^\ „„.uou.u • ( 3 - 8 ) 



k cosh(7rfc/2e) cosh(2A/c/e) 
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The form of In k L ' R (h) depends on the values of a and £. When either a < (L — 1)/2, 
or L/2 < a < (L + 2)/2 with £ < 0, ln^' R (n) is given by 

ln^> R (u) = -L [(oA + £)(tt-4£) + (|£| - 2A)vr + 2(2- a 2 )A 2 + 4u(A - u)] 



2vre 

cosh[(vr/2 - 2A)fc/e] ^ sinh[(vr/2 - 2X)k/e] 



-E 



oo 

ECOSIl[(J/( / Z — 4A)K/ fcj v 
fcsinh(7r/2e) ^ fccosh(7rfe/2e) 

cosh[(aA + £ - |f |)fc/e] cosh[(yr/2 - aX - f - |f |)jfe/e] cosh[(A - 2u)fc/e] 



| fc sinh(7r/c/2e) cosh(AA;/e) 

_ ^ cosh[(aA + g - \j\)k/e] sinh[(^/2 - aX - g - |f |)fc/e] cosh[(A - 2n)fc/e] 
^-^ fccosh(7rA;/2e) cosh(AA;/e) 

fc+i=l 

(3.9) 

On the other hand, when either a > (L + 3)/2, or L/2 < a < (L + 2)/2 with f > 0, 
hiK L ' R (u) is given by 

ln^. R ( n ) = J- [(oA+0(37r-4f) + (|£| -2A-vr)vr + 2(2-a 2 )A 2 + 4u(A-u)l 

^ cosh[(vr/2 -2X)k/e] ^ sinh[(vr/2 - 2X)k/s] 

4^ A;sinh(7rA;/2e) ^ k cosh(7r/c/2e) 

— fc+ 5 =i 

g cosh[(7r-aA-f-lf l)fc/e] cosh[(aA-7r/2+f- |f |)jfe/e] cosh[(A-2u)fe/e] 



k sinh(7rA;/2e) cosh(AA;/e) 
cosh\{n-aX-£-\£\)k/e\ sinh[(aA-7r/2+f- |f \)k/e] cosh[(A-2u)/c/e] 



fe=i 

oo 

y 

fecosh(7r/c/2e) cosh(A/c/e) 



(3.10) 



We note that once again Ihk s (u) is unchanged under the height reversal transfor- 
mation a — > L + 1 — a and f — > — f . 

The leading-order singular behaviour of the surface free energy is derived from 
the sums in equation (|3.8|). We find 



In k s 



ln(-i) L = 3 (mod 4) 
(-t) 71 "/^ otherwise. 



Hence in regime IV we have for all L 



vr L + 1 , 
2-« s = - = — . (3.12) 
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4 Band structure and ground state degeneracies 



In regimes III and IV the ABF models admit a number of coexisting phases [14, IT 



corresponding to ground state (zero temperature, \t\ — > 1) configurations on the 
lattice given by a pair (a, b) such that a — b = ±1. Explicitly, the ground states 
consist of a chequerboard arrangement with height a on one sublattice and height b 
on the other sublattice of the square lattice. At each boundary the boundary height 
a together with the choice of £ determines the favoured boundary configuration. 
If £ > the favoured boundary condition is {a, a + 1, a, a + 1, . . . }, and if £ < 
the favoured boundary condition is {a, a — 1, a, a — 1, ... }. Clearly, the boundary 
condition can pick out the phase near the boundaries. When the phase selected 
on the left and right boundaries is the same there will be a unique phase and a 
unique ground state. However, when the phase selected by the left boundary differs 
from that selected by the right boundary, there will be an interface formed between 
the two coexisting phases. In this case there will be an interfacial tension, the 
ground state will be highly degenerate corresponding to the many possible locations 
of the interface, and the double-row transfer matrix will have a continuous band 
of largest eigenvalues in the thermodynamic limit. More generally, we expect the 
spectrum of the double-row transfer matrix to consist of many such continuous 
bands of eigenvalues. These bands can overlap with each other but within each band 
the states are characterized by a fixed number of domain walls. Furthermore, the 
number of eigenvalues in each band can be determined by the simple combinatorics 
of counting the compatible ground state configurations. 

4.1 Regime III 

In regime III there are 2L — 2 coexisting phases. The corresponding ground state 
lattice configurations are given by the pairs (a, 6) with |a — b\ = 1, and all values 
1 < a, b < L allowed. Let us denote the left boundary condition with height a = a L 
by a + if £ L > and a~ if £ L < 0. Similarly, for the right boundary with height 
b = a R we have b + or b~ depending on the sign of £ R . This is convenient since under 
the restriction ( 1.15| ) the interfacial tensions and correlation lengths will depend only 



on the signs of £ L and £ R . The ground state degeneracies for a double-row transfer 
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matrix of N faces with fixed boundary heights a and b are then given by 



(4.1) 



+ ,+s - ,_v f(N+\a- b\)/2 

Af(a + ,b + ) = Af(a ,b )=( ^ & i 



For ease of reference, we define m(a, /3) and n(a, /3) by 

A^) ,4,) 

where a = a 1 * 1 and /3 = 6 =t . To see how these formulas arise, consider for simplicity 
an (a + , 6 _ ) boundary height configuration with b — a>2. Moving from a to b, there 
must be b — a — 1 ground state transitions (domain walls), starting with the state 
{a, a + 1, a, a + 1, . . . } and ending with the state {. . . , b— 1, b, b— 1, b}. If a and b are 
joined by a path of N steps, and we let mi,rri2, ■ ■ • ,m&_ _i be the steps at which 
the transitions occur, then the mi can take the values 

mi = 2,A,...,N + a-b + 2 (4.5) 
nii = mj_i + 1, mj_i +3, . . . , TV + a — 6 + i + l for 2 < i < b — a — 1. (4.6) 

We rewrite and sum over these combinations using the variables k±,k2, ■ ■ ■ fc(,_ a _i 
such that mj = 2/cj + i + 1. This gives 

jV + q-h N + a-b N+a-b 
fel=0 fe 2 =fcl fe(,_ a _i=fci,_ a _ 2 



in agreement with equation ( |4.2[ ). The other expressions may be derived from sym- 
metries under the interchange of a and b, and from the equivalence of (a~,/3) with 
iV steps and ((a - 1)+, (3) with TV + 1 steps. 

4.2 Regime IV 

First consider L even. In this case there are 2L — 4 phases, and the ground state 
configurations are as in regime III, with the exclusion of the (L/2,L/2 + 1) and 
(L/2 + l,L/2) pairs. In the case when L is odd, there are again 2L — 4 phases and 
the ground state lattice configurations are as in regime III, except now the pairs 
{(L - l)/2, (L + l)/2) and ((L + l)/2, (L + 3)/2) are excluded and replaced by the 
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disordered mixture ((L + l)/2, (L + l)/2 ± 1) where the heights on one sublattice 
can take the values (L — l)/2 and (L + 3)/2 independently at each site. 

For both even and odd L, we find that the number of eigenvalues in the first 
band is given as follows (|_xj denotes the largest integer less than or equal to x): 

• If 1 < a,b < [(L + l)/2\, excluding [(L + l)/2\ + , the degeneracies are the 
same as in regime III. 

• If [(L + 2)/2\ < a,b < L, excluding [(L + 2)/2\~, the degeneracies are the 
same as in regime III. 

• Otherwise, the degeneracies are 

AW)=A^n = ( (w ;^'- 2>/2 ) < 4 - 8 > 

5 Interfacial tensions 

The interfacial tension cr <a ' /3 ) between phases a = and (3 = b ± is given by 

_ a (a,P) = lim lim p-ilnTrA^)^ 2 (5.1) 

P^oo N^oo 



where 



A(«) = A(^(«) = r : ^'"^ , 1/2 (5.2) 



and a factor of the inverse temperature has been absorbed into the definition of 
<j( a 'P\ Here Dm^(u) denotes the maximum eigenvalue of the double-row transfer 
matrix D^ a ' a \u). Note that the bulk and surface contributions cancel out in the 
ratio A(u). From the inversion relation, crossing, and quasiperiodicity, it follows 
that the eigenvalues A(u) satisfy the simple functional equations 

A(u)A(u + A) = 1 A(u) = A(A - it) (5.3) 

subject to the periodicity 

e regime III 



A (it) = A(u + lire*) 



(5.4) 
2e regime IV. 
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In the thermodynamic limit the eigenvalues of A (it) form continuous bands. In 
this limit the sum in the trace of ( |5.1| ) is replaced by integrals over these bands. 
But in the limit of P large only the first band of largest eigenvalues is expected to 
contribute to the interfacial tension. From the counting of degeneracies, this first 
band contains N(a,(3) eigenvalues. The largest eigenvalue of D^ a '^(u) has zeros 
which accumulate on the lines Re(u) = —A/2 and Re(ti) = 3A/2 in both regimes, 
and also on the lines Re(it) = (3A — 7r)/2 and Re(u) = (it — A)/2 in regime IV. In 
addition, each eigenvalue in the first band of eigenvalues of D^ a '^\u) has n(a,f3) 
pairs of zeros in the strip —A/2 < Re(u) < 3A/2 at 

u = X/2±i(f) k k = 1,2, ... ,n(a,f3) (5.5) 

where each <p k is real and non-zero. 

We now seek a solution to the inversion relations ( [5.3| ) in the physical strip 
subject to the given zeros and periodicity. The required solution is 

n(o,/3) 

A(n)= ]J <S>(u + i(t> k )®{u-ict> k ) (5.6) 
k=i 



where 



and 



„Q (TTU 7T U|f\ 



(5.8) 



(L+l)/4 regime III 
(L+l)/2 regime IV. 

Recall that \t\ = \p 2 \ = exp(— 2tte) where the nome p is defined in equation (1.4). 

In the thermodynamic limit the distribution of each 4>k becomes dense, yielding 
a density p((j>i, ■ ■ ■ ,4> n (a,p))- So integrating over this band of eigenvalues for real 



values of u in the interval < u < A gives 

n(o,/3) 

Tr A(u) p / 2 ~ / ••• / p(d>i, d>„r n hi 



lim TrA(u) p/2 ~ / ••• / p((j>x, . . . , n(a m) TT \$>(u + i<f) k )\ p dfa . . . d<f) n{afi) . 
Jo Jo k=1 

(5.9) 



For P large this multiple integral can be evaluated by saddle point methods. The 
saddle point occurs at 

u + i<t> k = A/2 + i7re*/2 k = 1,2, ... ,n(a, (3). (5.10) 
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exp[-aW>] 



Hence 

' i? 4 (o,i*r) i" (a,/3) 

where k' is the conjugate elliptic modulus 

1 _p2n-l\ : 



fc'dtr)"^/ 2 (5.11) 



^ = 11(1+^1) • (5 - 12) 

n=l \ " / 

It follows that the corresponding critical exponent is given by 

/ m ((L + l)/4 regime III 

ct^'^ ~ |t|" n = v = { (5.13) 

({L+l)/2 regime IV. 

Notice that the phases are linearly ordered and that the only way to get from 
one phase to another is to pass through any intermediate phases. So let us define 
a (a,f3) _ f or non-ground state phases. Similarly, for L odd in regime IV, let 
us define er^ 1 ' -1 )/ 2 >( i + 1 )/ 2 ) = etc. for the disordered phase. Then the interfacial 
tensions are additive and it follows, for example, that if a < b, 

(a+,b+) - _((a+fc-l)+,(a+k)+) 



^ (T ((a+fc-l)+,(a+fc)+)_ ( 5M ) 



a 

k=l 

6 Correlation lengths 

Let <p(a) be a given function of the height a and let ao and ai be two heights 
in the same column separated by an even number £ of lattice spacings. Then, at 
large distances, the (truncated) pair correlation function decays exponentially with 
a correlation length £ given by 

(<p(a )(p(a£)) - (<p(a ))(<p(ae)) ~ exp(-^ _1 ). (6.1) 

By standard row transfer matrix arguments [GJ the correlation length can be calcu- 
lated from 

-C 1 = lim lim £~ 1 In V CjKAuf 12 (6.2) 

j 7^ max 



where the Aj(u) are eigenvalues of 



A(u) = —. — ^ (6.3) 
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and the c, denotes u- independent matrix elements. The sum excludes the maximal 
eigenvalue. Again the bulk and surface terms cancel out of the ratio A (it). From 
the inversion relation and crossing, it follows that the eigenvalues A(it) again satisfy 
the simple functional equations 

A(«)A(« + A) = 1 A(«)=A(A-«) (6.4) 



subject to the previous periodicity (|5.4|) . 

In the thermodynamic limit the sum in Q6.2| ) is dominated by the first band of 
(2 — 8 a i+ — S a l-)( N 2 2 ) eigenvalues. In this band the eigenvalues Aj(u) are analytic 
in an open strip containing < Re (it) < A, with two pairs of zeros at 

u = \/2±i<p 1 u = X/2±i(j)2. (6.5) 

The required solution of the inversion relation is therefore 

2 

A(u) = IJ$(u + ifo)$(u-»0 fc ) (6.6) 
fc=i 

where $(it) is given by fl5,7|) . Passing to the thermodynamic limit for real u in the 
interval < u < A we obtain 

lim V c^uf 2 - f 2 / 2 c(0i,02) JT + #i<#2- (6.7) 

Carrying out the saddle point analysis as before gives 

exp(-r 1 ) = fc'(IC) (6-8) 
where k' is the conjugate elliptic modulus and hence 

t |( L + 1 )/ 4 regime III 

4 ~ |t| z/ = < (6.9) 

l(JD + l)/2 regime IV. 

It can be seen that, after allowing for differences in formulation, all of the above 
results for the interfacial tensions and correlation lengths agree with the known 
results in the case of the Ising model Q (L = 3) and the interacting hard square 
model [|l6| (L = 4). Notice also that we have the simple relation 

r i = 2a( 1+ < 2+ ). (6.10) 

The factor of 2 is easily understood because in the saddle point integrals the con- 
tribution to the inverse correlation length involves two domain walls, whereas the 
contribution to this interfacial tension derives from just one such domain wall. 



17 



Finally, we note that in regime III the critical exponents a, a s , fx and v satisfy 
the scaling relations |], [l?], [l8| 



a 



2u 



a + ia 



(6.11) 



However, these scaling relations break down in regime IV, since the regime IV line 
of exact solution approaches the tricritical point t = at a tangent to the critical 
line PI. 



7 Discussion 



In this paper we have calculated the surface free energies, interfacial tensions and 
correlation lengths of the ABF models in regimes III and IV, all by solving a rel- 
atively simple inversion relation. The methods are general and can be applied to 
other solvable lattice models such as the CSOS models and dilute lattice mod- 
els. The methods employed here could be applied to regimes I and II with fixed 
boundary conditions where the heights alternate along the boundary. However, in 
regimes I and II, the coexisting phases are selected out by a saw-tooth variation in 
the ground state configuration [14, and unfortunately these saw-tooth variations 



on the boundary cannot be handled within our formalism at present. So, strictly 
speaking, the current methods are best suited to systems that exhibit at most two 
independent sublattices in the structure of their ground states. 

Of course, the correlation length along the strip is a bulk property and not 
a surface property and so its calculation should not necessitate the introduction 
of a boundary. Accordingly, the generalized inversion relation method used here 
to obtain the correlation length can be applied in the case of periodic boundary 
conditions simply by working with a double-row transfer matrix 

D{u) = T(u)T{X-u). (7.1) 

Since the correlation length is independent of the spectral parameter u, and since 
D(u) and T(u) yield the same correlation length at the isotropic point u = A/2, the 
two transfer matrices must lead to the same correlation length for all values of u. 
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